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CN Abstract 

Bernstein-Sato polynomial of a hypersurface is an important object with numerous applica- 
^ tions. It is known, that it is complicated to obtain it computationally, as a number of open 

^ questions and challenges indicate. In this paper we propose a family of algorithms called 

^ checkRoot for optimized check of whether a given rational number is a root of Bernstein- 

p—i Sato polynomial and the computations of its multiplicity. This algorithms are used in the new 

approach to compute the whole global or local Bcrnstein-Sato polynomial and 6-function of 
a holonomic ideal with respect to weights. They are applied in numerous situations, where 
there is a possibility to compute an upper bound for the polynomial. Namely, it can be 
achieved by means of embedded resolution, for topologically equivalent singularities or using 
^ the formula of A'Campo and spectral mimbers. We also present approaches to the logarith- 

mic comparison problem and the intersection homology D-modulc. Several applications are 
^ presented as well as solutions to some challenges which were intractable with the classical 

^ methods. One of the main applications consists of computing of a stratification of affine 

space with the local 6-function being constant on each stratum. Notably, the algorithm we 
propose docs not employ primary decomposition. Also we apply our results for the compu- 
tation of Bernstein-Sato polynomials for varieties. The methods from this paper have been 
^ implemented in Singular:Plural as libraries dmod.lib and bfun.lib. All the examples 

O from the paper have been computed with this implementation. 
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1. Introduction 



Through the article we assume K to be a field of characteristic 0. By i?„ we denote the 
ring of polynomials K[xi, . . . in n variables over K and by Dn we denote the ring of 
K-linear partial differential operators with coefficients in i?„, that is the n-th Weyl algebra 
[19] . The ring Dn is the associative K-algebra generated by the partial differential operators 
di and the multiplication operators Xi subject to relations 

{diXj = Xjdi + 6ij,XjXi = XiXj, djdi = didj \ I < i, j < n}. 

That is, the only non-commuting pairs of variables are they satisfy the relation 

idi + 1. We use the Lie bracket notation [a, b] := ab — ba for operators a, 6, then e.g. 
the latter relation can be written as = 1. 

Finally, we denote by Dn[s\ the ring of polynomials in one variable s with coefficients in 
the n-th Weyl algebra, i.e. Djys] = Dn ®k 

Let us recall Bernstein's construction. Given a non-zero polynomial / G Rn in n variables, 
we consider M = Rn[s, j] ■ which is by definition the free Rn[s, j]-module of rank one 
generated by the formal symbol /'^. Then M has a natural structure of left -module. 
Here the differential operators act in a natural way, 

9.(^(.,x)-r) = (^ + .^(.,x)^^)-r eM (1) 

Theorem 1.1 (Bernstein |3]). For every polynomial f E Rn there exists a non-constant 
polynomial b{s) G K[s] and a differential operator P{s) G -D„[s] such that 

P{s)f-r = b{s)-f G Rn[s,j]-f = M. (2) 

The monic polynomial b{s) of minimal degree, satisfying ^ is called the Bernstein-Sato 
polynomial or the global 6-function. 

This paper is organized as follows. In Section |2| the checkRoot family of algorithms 
for checking rational roots of the global and local Bernstein-Sato polynomial is developed. 
We also show how to compute the 6-function of a holonomic ideal with respect to a certain 



weight vector. In Section 3.1 we show how to obtain an upper bound in various situations (by 
using an embedded resolution, for topologically equivalent singularities, by using A'Campo's 
formula and spectral numbers). In particular, we demonstrate a complicated example of 
(non-isolated) quasi-ordinary singularity. 

In Section |4.2 we discuss the possibilities to obtain integral roots of the 6-function and 



apply it to the computation of the minimal integral root in the context of Intersection Ho- 
mology D-module and Logarithmic Comparison Theorem. In Section [5] we present a new 
method for computing the stratification of affine space, according to local Bernstein-Sato 
polynomials. 

We want to stress, that Bernstein-Sato polynomials for most of the examples, presented 
in this paper, cannot be computed by direct methods with any computer algebra system 
including SINGULAR: PLURAL [TU]. Indeed, these examples were known as open challenges 
in the community and here we present their solutions for the first time. 
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The timings for examples in this paper were performed on a PC with 4 Dual Core AMD 
Opteron 64 Processor 8220 (2800 MHz) (only one processor available during the computation) 
equipped with 32 GB RAM (at most 16 GB available during the computation) running 
openSUSE 11 Linux. 

2. The checkRoot Family of Algorithms 

For the sake of completeness, some of the ideas coming from [H], as well as some results 
and their proofs have been included here. 

Several algorithms for computing the 6-function associated with a polynomial are known, 
see e.g. |2SII2S1EZ], |2H], [32], [II]- However, from the computational point of view 

it is very hard to obtain this polynomial in general. Despite significant recent progress, only 
restricted number of examples can be actually treated. In order to enhance the computation 
of the Bernstein-Sato polynomial via Grobner bases, we study the following computational 
problems. 

1. Obtain an upper bound for 6/(s), that is, find B{s) G K[s] such that &/(s) divides B{s). 

d 

B{s) = His - 
1=1 

2. Check whether is a root of the 6-function. 

3. Compute the multiplicity of as a root of bf{s). 

There exist some well-known methods to obtain an upper bound for the Bernstein-Sato 
polynomial of a hypersurface singularity once we know, for instance, an embedded resolution 
of such singularity [12]. However, as far as we know, there is no algorithm for computing the 
6-function from this upper bound. In this section we present algorithms for checking whether 
a given rational number is a root of the 6-function and for computing its multiplicity. As 
a first application, using this idea, we could obtain 6/(s) for some interesting non-isolated 



singularities, see Example 3.3 below. 



From the definition of the 6-function it is clear that 

(6K^)) = (Annz,„w(r) + (/))nK[s]. (3) 

In fact, this is another way of defining the Bernstein-Sato polynomial. This equation was 
used to prove the main result of this section, namely Theorem 2.1[ 



Theorem 2.1. Let R be a IK-algebra, whose center contains K[s]. Let q{s) G K.[s] be a 
polynomial in one variable and I a left ideal in R satisfying / fl K[s] ^ 0. The following 
equalities hold: 

1. (J + R{q{s))) n K[s] = / n K[s] + K[s] {q{s)), 

2. (/ : q{s)) n K[s] = {10 K[s]) : q{s), 

3. (/ : g(s)°°) n K[s] = (/ H K[s]) : g(s)°°. 

In particular, using I = Ann£)^[s]if^) + (/) ^ -D„[s] in the above equation we have 
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[Annn„i4n + Dn[s]{f,q{s))] nK[s] = {bf{s),q{s)) = { gcd{bf{s),q{s))), 

gcd{bf{s),q{s)) 



[(Annz,„w(r) + D^Mf)) ■ ?(^)] H = : q{s) = { -.r^fl^. >, 



. [(Annz,„w(r) + D4s]{f)) : n K[s] = (bfis)) : g(.)-. 

Proof. Let 7^ be a generator of /nK[s]. At first, suppose that h{s) E (/+i?(g(s)))nK[s]. 
Then we have 

h{s)=Pis) + Qis)qis) (4) 

where P{s) G / and Q{s) G -R. Let (i(s) be the greatest common divisor of b{s) and q{s). 
There exist and qi{s) such that d{s)bi{s) = b{s) and d{s)qi{s) = q{s), and hence 

bi{s)q{s) = qi{s)b{s). Since s commutes with all elements in R, multiplying the equation Q 
by bi{s), one obtains 

bi{s)h{s) = bi{s)P{s) + Q{s)qi{s)b{s) G / 

Thus, bi{s)h{s) G /nK[s] = (6(s)) and therefore h{s) G : = (ci(s)) = /nK[s] + 

(g(s)). The other inclusion follows obviously. The second and the third parts can be shown 
directly and the proof is complete. □ 

Note that the second (resp. third) part of the previous theorem can be used to heuris- 
tically find an upper bound for bf{s) (resp. the roots of bf{s)). Since q{s) is in the center 
of -D„[s], the quotient and saturation ideals can be computed effectively via the kernel of a 
module homomorphism procedures, cf. [T3]. More classical but less effective approach is to 
use the extra commutative variable, say T, and the formula 

/ : q{sr = D^[s, T]{1, 1 - Tq{s)) H D^[s]. 



Let us see an example to illustrate how useful could be Theorem 2.1 



Example 2.2. Let / G C[x,|/] be the polynomial x{x'^ + y^). The annihilator of/* in D[s] can 
be generated by the operators -Pi(s) = ?>xy'^dx — y^dy — 3x'^dy and P2{s) = Sxd^: + "iydy — 9s. 
Consider the univariate polynomial 

q{s) = {s + l){s + 5/9)(s + 8/9)(s + 10/9)(s + 7/9){s + ll/9)(s + 13/9). 

Computing a Grobner basis, one can see that the ideal in D[s, T] generated by {Pi{s), P2{s), 
/, 1 — Tq{s)} is the whole ring. From Theorem 2.1 (3), one deduces that q{s) contains all the 
roots of bf{s). Using this approach we only have to check whether an ideal is the whole ring 
or not. Therefore any admissible monomial ordering can be chosen, hence the one, which is 
generically fast. 

Given an arbitrary rational number a, let us consider the ideal la C Dn[s] generated by 
the annihilator of the polynomial / and s + a. Theorem 2.1 (1) says that the equality 
la = Dn[s] holds generically (this is clarified in Corollary 2.3 below). Hence the roots of the 
Bernstein-Sato polynomial are the rational numbers for which the condition Jq 7^ Dn[s] is 
satisfied. This allows one to work out with parameters, that is over ^{a){x,dx)[s], and find 
the corresponding complete set of special parameters. The latter procedure is algorithmic [15] 
and implemented in SINGULAR. Note, that the set of candidates to obstructions, returned 
by the latter algorithm is in general bigger, than the set of real obstructions. 
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Corollary 2.3. Let {Pi{s), . . . ,Pk{s)} be a system of generators of the annihilator of in 
Dn[s\. The following conditions are equivalent: 

1. a G Q>o is a root ofbf{—s). 

3. Dn{Pii-a),...,Pk{-a),f) ^ Dn. 
Moreover, in such a case Dn[s]{Pi{s), . . . , -Pfc(s), /, s + a) fl K[s] = K[s]{s + a). 

Proof Take J = D„[s](Pi(s),...,Pfc(s),/,s + «) and K = JnD„ = (Pi(-«), . . . Pfc(-a), /). 
Since 

J = D[s]^ Jn K[s] = K[s] ^ K = Dn, 
and gcd(6/(s), s + a) = 1 if and only if bf{—a) ^ 0, the result follows from applying Theorem 



2.1 using = s + a. □ 



Once we know a system of generators of the annihilator of in -D„[s], the last corollary 
provides an algorithm for checking whether a given rational number is a root of the 6-function 
of /, using Grobner bases in the Weyl algebra. 

Algorithm 1 CHECKRootI (checks whether a G Q>o is a root of bf{—s)) 



Input 1 
Input 2 
Output 



{Pi(s), . . . , Pfc(s)} C P)„[s], a system of generators of Ann£,„[^](/" 
/, a polynomial in P„; a, a number in Q>o; 
true, if a is a root of bf{—s); false, otherwise; 



n 



K := (Pi(-a), . . . , Pfc(-«), /); >K = JnD^CD 

G := reduced Grobner basis of K w.r.t. ANY term ordering; 
return {G ^ {1}); 

2.1. Multiplicities 

Two approaches to deal with multiplicities are presented. We start with a natural gener- 



alization of Corollary 2.3 



Corollary 2.4. Let rria be the multiplicity of a as a root of bf{—s) and let us consider the 
ideals Ji = Ann£)^^[5](/*) + (/, (s + ay~^^) C i = 0, . . . ,n. The following conditions are 

equivalent: 

1. rria > i. 

2. Ji, nK[s] = {{s + ay+^). 

3. (s + aY i Ji. 

Moreover if Dn\s\ 2 2 2 ' ' ' 2 Jm-\ = Jm, then rua = m. In particular, m < n and 

Jm~l Jm ' ' ' Jfi' 

Proof. 1 <^==^ 2. Since rUa > i ii and only if gcd(6/(s), (s + a)*"*"^) = {s + aY^^, the equivalence 



follows by applying Theorem 2.1 (1) using q{s) = (s + ay^^ 



2^3. If (s + a)' e JiH K[s], then clearly Ji n K[s] 2 {{s + 
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3 =^ 2. Let h{s) G K[s] be the monic generator of the ideal Jj fl K[s]. Since (s + a)*"*"^ G 
Jj n K[s] = {h{s)), there exists j < i + 1 such that h{s) = (s + a)-'. Suppose that j < i. 
Then, 

(s + ay = {s + ay-^{s + ay = {s + ay-^h{s) E J,. 

That, however, contradicts 3 and thus j = i + 1. 

The rest of the assertion follows by applying the above result using i = m and i = m — 1, 
since (s + a)™ G Jm and (s + a)"^^^ ^ Jm-i from the hypothesis. □ 

Again once we know a system of generators of the annihilator of in -D„[s], the last 
corollary provides an algorithm for checking whether a given rational number is a root of 
the 6-function of / and for computing its multiplicity, using Grobner bases for differential 
operators. 



Algorithm 2 checkRoot2 (computes the multiplicity of a G Q>o as a root of s)) 
Input 1: {Pi(s), . . . , -Pfc(s)} C Z}„[s], a system of generators of Ann£)^[s](/*); 
Input 2: /, a polynomial in a, a number in Q>o; 

Output: rria, the multiplicity of a as a root of s); 

for i = to n do 

J:=D^[s]-{P^is),...,P,is)J,is + ay+'y 
G := Grobner basis of J w.r.t. ANY term ordering; 
r := normal form of (s + a)* with respect to G; 
if r = then 
rria '■= i] 
break 
end if 
end for 
return m^; 



t> r = ^ {s + ay e Ji 
> leave the for block 



Proof, (of Algorithm [2]) . 

Termination: The algorithm checkRoot2 clearly terminates and one only has to consider 
the loop from to n because the multiplicity of a root of &/(s) is at most n, see [30] . 
Correctness: Corollary 2. 4| implies the correctness of the method. □ 

Remark 2.5. There exists another version of checkRoot2 with just one step, due to the 
formula, see Corollary above, 

{Annn^[s]{n + Dn[s]{f,{s + ar))nK[s] = {{s + ar-). 

However, this method only seems to be useful when the multiplicity is close to n, otherwise 
checkRoot2 is more effective. The reason is that in general, the multiplicity is far lower than 
the number of variables. 

This algorithm is much faster, than the computation of the whole Bernstein polynomial 
via Grobner bases, because no elimination ordering is needed for computing a Grobner basis 
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of J . Also, the element (s + a)*"*"^, added as a generator, seems to simplify tremendously such 
a computation. Actually, when i = it is possible to eliminate the variable s in advance and 



we can perform the whole computation in see Corollary 2.3 (3) above. 

Nevertheless, Algorithm [2] meets the problem to calculate on each step a Grobner basis 
Gi for an ideal of the form / + ((s + aY^'^) and the set Gj-i is not used at all for such 
computation. A completely new Grobner basis has to be performed instead. The classical 
idea of quotient and saturation are used to solve this obstruction. In particular, the following 
result holds. 

Corollary 2.6. Let rria be the multiplicity of a as a root of bf{—s) and let us consider the 
ideal I = Ann [,,](/'') + Dn[s]{f) . The following conditions are equivalent: 

1. nia > i- 

2. {I : (s + aY) +D4s]{s + a) ^ D^s]. 

3. {I ■.{s + ay)\s=-a^Dn. 

Proof. Given J C Dn[s] an ideal, we denote by bj{s) the monic generator of the ideal Jn]K[s]. 
Then, from Theorem 2.1 (1), condition 2 is satisfied if and only if —a is a root of &/:(s-j-av(s). 



This univariate polynomial is nothing but bf{s)/ gcd(6/(s), (s + a)*), due to Theorem 2.1 (2). 
Now the claim follows from the obvious equivalence 



nia > i <^=^ (s + a) 



bfis) 



gcd(6/(s),(s + a)0' 

□ 



Since s+a belongs to the center of -D„[s], the ideal / : (s+a)' can recursively be computed 
by the formulas 

I: (s + a) = {lr]Dn[s]{s + a))/{s + a), 
I:{s + ay = {I : {s + ay~^) : (s + a). 

The following is a sketch of another algorithm for computing multiplicities using quotient 
ideals. The termination and correctness follow from the above corollary. 

Algorithm 3 checkRootS (computes the multiplicity of a G Q>o as a root of &/(— s)) 
Input 1: {Pi(s), . . . , Pk{s)} C a system of generators of Ann£i^[s](/*); 

Input 2: /, a polynomial in a, a number in Q>o; 

Output: rua, the multiplicity of a as a root of s); 

m:=0; / := (Pi(s), . . . , Pfc(s), /); J := / + (s + «); 
while G ^ {1} do 
m := m + 1; 

I ■= I : {s + a); t> I : {s + af 

J := I + Dr,[s]{s + a); (or J := 

G := reduced Grobner basis of J w.r.t. ANY term ordering; 
end while 
return m; 
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Remark 2.7. Several obvious modifications of tlie presented algoritlims can be useful depend- 
ing on tlie context. Assume, for instance, that q{s) is a known factor of the Bernstein-Sato 
polynomial and one is interested in computing the rest of bf{s). Then the ideal / : q{s) 
contains such information. This easy observation can help us in some special situations. 

Remark 2.8. Define the reduced Bernstein-Sato polynomial of / G -R„ to be = 

bfi^)/i^ + !)• Recall, that the Jacobian ideal of / is J/ = (^, . . . , ^) C K[x]. It is 
known, that taking (/) + Jf instead of (/) has the following consequence 

(Ann^W r + (/, 1^, • • • , S:)) n ns] = {b'jis)) = (^). 

Hence, all the algorithms above can be modified to this setting, resulting in more effective 
computations. This is the way it should be done in the implementation. We decided, however, 
not to modify the description of algorithms in order to keep the exposition easier. 

2.2. Local versus global b-functions 

Here we are interested in what kind of information one can obtain from the global b- 
function for computing the local ones and conversely. In order to avoid theoretical problems 
we will assume in this paragraph that the ground field is C. 

Several algorithms to obtain the local 6-function of a hypersurface / have been known 
without any Grobner bases computation but under strong conditions on /. For instance, 
it was shown in [17] that the minimal polynomial of —dtt acting on some vector space of 
finite dimension coincides with the reduced local Bernstein polynomial, assuming that the 
singularity is isolated. 

Remark 2.9. Recall, that the singular locus of V{f) is V{{f, . . . , ^)). One can define 
the local 6-function or local Bernstein-Sato polynomial as follows. Let p G be a 
point and trip = {{xi — pi, . . . ,Xn — Pn}) C Rn the corresponding maximal ideal. Let Dp be 
the local Weyl algebra at p, that is Weyl algebra with coefficients from C[xi, . . . , Xn]p instead 
of Rn = C[xi, . . . ,Xn]- From the Bernstein's functional equation ^ it follows that 3-P(s) G 
D[s],bf{s) G IK[s], such that P{s)f ■ = b{s) ■ holds. Hence, since over C[xi, . . . ,Xn]p 
we have got polynomial invertible, there exist 3-Pp(s) G Dp[s],bf^p{s) G K[s], such that 
Pp{s)f ■ = bf^p{s) holds. We define local Bernstein-Sato polynomial to be the univariate 
monic polynomial &/,p(s) of the minimal degree, such that the above identity holds. 

Theorem 2.10. (Briangon-Maisonobe (unpublished), Mebkhout-Narvdez fW^) Let fe/,p(s) 
the local b- function of f at the point p G C" and bf{s) the global one. Then bf{s) = 
Icmpec" bf^p{s) = lcmpgs(/) bf,p{s). 

The previous Theorem can be very useful for computing the global 6-function using the 
local ones. Let us see an example. 

Example 2.11. Let C the curve in given by the equation / = (y^ — x^)(3x — 2?/ — l){x + 2y). 
This curve has three isolated singular points (0,0), (1, 1) and (1/4, —1/8). The following is 
its real picture. 
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Figure 1: The cup (2,3) with two hues. 



The library gmssing.lib contains a procedure bernstein, which computes the local 
6-function at the origin. Moving to the corresponding points we can also compute &/,pj(s). 

bf,p^{s) = (s + l)2(s + 5/8)(s + 7/8)(s + 9/8)(s + ll/8) 
bf,P2is) = (3 + l)^(s + 3/4)(. + 5/4) 
hpM = is + l)'is + 2/3)is + A/3) 



From this information and using Theorem 2.10 the global 6-function is 



{s + l)\s + 2/3)(s + 5/8)(s + 3/4)(s + 7/8)(s + 4/3)(s + 5/4)(s + 9/8)(s + 11/8) 



The computation of the global 6-function with Theorem 2.10 is effective, when the singular 
locus consists of finitely many isolated singular points. The Singular library gmssing.lib 
implemented by M. Schulze [31] and based on his work [32] allows one to compute invariants 
related to the the Gauss-Manin system of an isolated hypersurface singularity. In the non- 
isolated case the situation is more complicated, since no Gauss-Manin connection exists. For 
computing the local 6-function in this case (which is important on its own) we suggest using 
the global 6-function as an upper bound and a local version of the checkRoot algorithm, see 
Section 12.2.21 below. 

In [21], H. Nakayama presented an algorithm for computing local 6-functions. One step 
in his algorithm uses a bound for the multiplicity of a given rational root of the global b- 
function. Then the algorithm checks if this multiplicity agrees with the local one. This 
approach is very similar to our checkRoot algorithm. 

2.2.1. Localization of non- commutative rings 

We recall some properties of rings of fractions in non-commutative setting. The reader is 
referred to |9| and [T9] for further details. 



Definition 2.12. Let i? be a ring and S ^ R a multiplicatively closed set. A left ring of 
fractions for R with respect to S" is a ring homomorphism cf) : R ^ Q such that: 

1. 0(s) is a unit of Q for all s G S*. 

2. Each element of Q has the form for 0(s)^^0(r) for some r & R and s G S". 

3. ker(0) = {r ^ R \ sr = for some s G S}. 
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Right rings of fractions are defined analogously. 



Theorem 2.13. There exists a left ring of fractions for R with respect to S if and only if S 
is a left denominator set, that is, the following conditions hold: 

• Left Ore condition: for each r E R and s E S, there exist r' E R and s' E S such that 
s'r = r's, that is, Sr H Rs 

• Left reversible: if rs = for some r E R and s E S, then 3r' G R such that sr' = 0. 

In such a case, the pair {Q, (p) is universal for homomorphisms if : R T such that ^p{S) 
consists of units of T and therefore Q is unique up to unique isomorphism. Moreover, if R 
also has a right ring of fractions Q' with respect to S then Q Q' . 

Because of the uniqueness, the left ring of fractions Q (when it exists) is often denoted 
by S^^R, and the natural map (j) : R ^ S^^R is called the localization map. To simplify 
notation the elements of S~^R are denoted by s~^r, even when ker0 ^ 0. Two quotients 
s^^ri and S2^r2 are equal if and only if there exist s E S and a E R such that asi = SS2 and 
ari = sr2. The localization for left (resp. right) modules can be generalized in the obvious 
way and it is verified S-^M = S-^R 0r M (resp. MS-^ = M ®r RS-^). 

Recall the following two classical results on localizations. 

Lemma 2.14. Let Ri ^ R2 he a ring extension and S G Ri a multiplicatively closed set. 
Assume S~^Ri and S~^R2 exist and consider the corresponding localization maps 0i : i?i — )■ 
S~^Ri and 02 : -R2 — ^ S~^R2. Let j : S~^Ri — i- S~^R2 be the map induced by i. Then j is 
infective and for every left ideal I ^ R2 one has 

s-'^ins-^Ri = s-\inRi). 

Note that in the previous lemma, 5*"^/ is the extension of / to S~^R2 while S~^{I fl -Ri) 
is the extension of / fl -Ri to S~^Ri. 

Lemma 2.15. Let R be a ring, S G R a multiplicatively closed set and I G R a left ideal. 
Assume S'^R exists. Then S~^I is not the whole ring S~^R if and only if I (1 S = ^. 

Example 2.16. Let R = D he the classical ra-Weyl algebra and 5* = K[x] \ m.p, where p E K" 



is an arbitrary point, cf. Remark 2.9 Then S* is a left and right denominator set as in the 



statement of Theorem 2.13, and the localization (K[x] \ trip) is naturally isomorphic to 



Dp. Analogous construction also holds for the extension D[s] = K[s] ®k -D. 
2.2.2. Local version of the checkRootl algorithm 



Theorem 2.1 is general enough to apply also for checking rational roots of local Bernstein- 
Sato polynomials. To simplify the exposition, we concentrate our attention on the local 
version of checkRootl algorithm. See Section |5] for other generalizations. 

Let / G C[xi, . . . , x„] be a polynomial, p G C" and a G Q. Then the first part of Theorem 



2.1, see also Corollary 2.3, tells us that (s + a) is a factor of the local 6-function at p if and 
only if the left ideal 

AnnB^[,](r) + DpM(/,s + «) (5) 
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is not the whole ring -Dp[,s]. From Lemma 2.14 using Ri = D[s], R2 = D{t,dt) := D 



K{t, dt I dt ■ t = t ■ dt + 1), S = C[x] \ trip and / = Anno(j,at)(/'') = // the Malgrange ideal 
associated with /, one obtains 

AnnD,[s]{r) = Dp[s] Annz5w(/'). 

Proposition 2.17. Let {Pi{s), . . . ,Pk{s)} be a system of generator 0/ Ann£)[s](/'') and con- 
sider the ideal I = D[s] (-Pi(s), . . . , Pk{s) , f , s + a) . Then we have 

{s + a)\bf^p{s) ^ pGF(/nC[x]). 

Proof. From the above discussion, Dp[s]I equals the ideal of the equation ([s]) and thus (s + a) 
is a factor of iff Dp[s]I 7^ -Dp[s]. Now, by Lemma 2.15 using R = D[s] and S = C[x]\mp 



Dp[s]I Dp[s] ^ /n (C[x] \mp) = ^ /nC[x]Cmp 

and the claim follows. □ 

There are several ways to check whether an ideal / C Dp[s] is proper or not. However, 
it is an open problem to decide which one is more efficient. Mora division and standard 
bases techniques seem to be more suitable in this case, since otherwise a (global) elimination 
ordering is needed. On the other hand, using this approach, such orderings are unavoidable 
for obtaining the stratification associated with local 6-functions, see Section [5] where several 
examples are showed. 

2.3. h-functions with respect to weights and checkRoot 

The 6-function associated with a holonomic ideal with respect to a weight is presented. 
We refer [28] for the details. Let 7^ G ]R>o and consider the \^-filtration with respect 
to w, {Vm \ m E Ij} = V on D where is spanned by | —wa + wf3 < m} over K. 

That is, Xi and di get weights —Wi and Wi respectively. Note, that with respect to such 
weights the relation diXi = Xidi + 1 is homogeneous of degree 0. The associated graded ring 
gr^(D) = Kra/Kri-i is isomorphic to D, which allows us to identify them. 

For a non-zero operator 

P= a^fsx^'d'^ e D, 

the maximum maxa^/jj— + w/3 \ 7^ 0} G M is denoted by ord^(P) and the principal 
symbol of P is the V^-homogeneous operator given by 

a^(P) := (^c.px'^d^. 

— uiQ+M)/3=ordv (P) 

Additionally, for a given ideal I ^ D, the associated graded ideal is defined as the vector 
space spanned by all its principal symbols, that is, gr^(/) := IK ■ {cr^(P) | P G /}. 

Sometimes, the principal symbol (resp. associated graded ideal) is called the initial form 
(resp. initial ideal) and it is denoted by in(_u,^^)(P) (resp. in(_^_^)(J)). 
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Definition 2.18. Let / C -D be a holonomic ideal. Consider 7^ w G 



J2i=i "^i^idi- Tlien gr^(/) nK[s] 7^ is a principal ideal in 
the global b-function of I with respect to the weight w. 



^>o 



and s 



[s] . Its monic generator is called 



Although Theorem 2.1 can not be applied in this setting, since s = ^^WiXidi does not 



belong to the center of the algebra, a similar result still holds, due to the properties of the 



filtration, see Proposition 2.19 below. Also Corollaries 2.3, 2.4 and 2.6 can be established 
using initial parts instead of annihilators. 

Proposition 2.19. ( gr^(/) + gr^(L))(g(s)) ) n K[s] = gr^(/) n K[s] + K[s](g(s)). 

Proof. Actually it is an easy consequence of being treated with ^^-homogeneous ideals. Con- 
sider h{s) = Q + R ■ q{s), where Q G gr^(J) and R G gr^(Z)). Taking ^^-homogeneous 
parts in the above expression, one finds Qq G gr^(/) and Rq G gr^(D) of degree such that 



his) 



Theorem 2.1 



Qo + Rq- q{s). Now, since q{s) commutes with Qq, one can proceed as in the proof of 
!)• 



□ 



Many algorithms in the realm of D-modules are based on the computation of such b- 
functions. For some applications like integration and restriction, only the maximal and the 
minimal integral roots have to be computed. 

However the above proposition can not be used to find the set of all integral roots, since 
neither upper nor lower bound is known in advance. For instance, N. Takayama used the 
following easy example to show the general unboundness: / = {xdi + k), k G Z is Di- 
holonomic and in(_i^i)(/) fl C[s] = {s + k) with s = tdt- 



3. Computing b-functions via Upper Bounds 

As different possible ways to find upper bounds, we present embedded resolutions, topo- 
logically equivalent singularities and A'Campo's formula. Depending on the context local or 
global version of our algorithm is used. 

3.1. Embedded resolutions 

In this part of the paper we will work over the field C of the complex numbers. However, 
in actual computation we can assume that the ground field is generated by a finite number 
of (algebraic or transcendental) elements over the field Q of the rational numbers and that 
the algebraic relations among these elements are specified. 

Definition 3.1. Let /i : F — )■ C" be a proper birational morphism. We say that /i is a global 
embedded resolution of the hypersurface defined by a polynomial / G C[x], X = V{f), if the 
following conditions are satisfied: 

1. y is a non-singular variety. 

2. h : Y \ h~^{X) — C" \ X is an isomorphism. 

3. h~^{X) is a normal crossing divisor. 
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Since h~^{X) is a normal crossing divisor, the morphism F = foh: Y^Cis locally 
given by a monomial. Hence, we can define the 6-function of F as the least common multiple 
of the local ones. If F is locally given by the monomial at the point p, then 

i=l ^ j=l ^ l<«j<aj \<k<n ^ 

The following is the global version of the classical result by Kashiwara The upper 
bound statement is due to Varchenko ([SB]) and Saito ([SHIED])- 

Theorem 3.2. For f G Rn, there exists an integer k such that bf{s) is a divisor of the 
product bp{s)bF{s + 1) ■ ■ ■ bpis + k). Moreover < k < n — 1. 

Proof. Since /i is a global embedded resolution of X = V{f), h induces a local embedded 
resolution of the germ {X,p) at every point p G X. Now, the existence of k > with the 
divisibility property follows from the theorem by Kashiwara [12j and from the fact that the 



global 6-function is the least common multiple of the local ones, see Theorem 2.10 The proof 



for the upper bound can be found in the references above. □ 

This theorem allows one to find upper bounds also for the global case. Let us see an 
example to show how one can apply the algorithm checkRoot in order to compute the b- 
function. 

Example 3.3. Let / = {xz + y){x'^ + + xy^) G Q[x,?/, z] and -Bi(s) = 6^.5(3)6^,18(5)6^24(5). 



Since every root of bf{—s) belongs to the real interval (0,3), see Theorem 3.2, computing 
an embedded resolution of the singularity and using Kashiwara's result [12], we obtain that 
B{s) = Bi{s)Bi{s + l)Bi{s + 2) is an upper bound for 6/(s). Once we know a system of 
generators of the Ann£)^[s] checking whether each root of the upper bound is a root of the 
Bernstein-Sato polynomial was easy. It took less than 5 seconds except for those ones which 
appear in the table below. We also observe that when a candidate is not a root indeed, the 
computation is very fast. To the best of our knowledge, this example (first appeared in [8]) 
is intractable by any computer algebra system. 

bj{s) = (s + lf{s + 17/24)(s + 5/4)(s + ll/24)(s + 5/8)(s + 31/24)(s + 13/24) 
(s + 13/12)(s + 7/12)(s + 23/24)(s + 5/12)(s + 3/8)(s + ll/12)(s + 9/8) 
(s + 7/8) (s + 19/24) (s + 3/4) (s + 29/24) (s + 25/24) 

The running time is given in the format minutes: seconds. 



Root of 5(-s) 


Running time 


Root of 6/(-s) ? 


5/4 


29:16 


Yes 


31/24 


26:16 


Yes 


29/24 


7:51 


Yes 


9/8 


0:35 


Yes 
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TT 



Let us give a brief indication for computing a global embedded resolution of /. Consider 
C'^ — )■ C^, the blow-up of with center in Z = {x = y = 0}. Denote Vi = V{xz + y) 
and V2 = V{x^ + y^ + xy'^) the two components of V{f) and Vi, V2 their corresponding 
strict transforms. The exceptional divisor E has multiplicity 5 and Vi and V2 do not meet. 
Moreover Vi and E intersect transversally. The local equation of V2 U is given by the 
polynomial y^{x'^ + y + xy). Now, one can proceed as in the case of plane curves, since 
the local equation involves just two variables. Finally, we obtain seven divisors with normal 



crossings, see Figure 3.1 , This method can also be applied to the family {xz + y)g{x, y) under 
some extra conditions on g{x,y). 




Figure 2: Embedded resolution of V{{xz + y)(x'^ + + xy'^)) 



Remark 3.4. To the best of our knowledge, resolution of singularities has never been used 
before for computing Bernstein-Sato polynomials in an algorithmic way. Recall that an 
embedded resolution can be computed algorithmically in any dimension and for any affine 
algebraic variety [5]. 

One can find upper bounds for the case of hyperplane arrangements by computing an 
embedded resolution. This allows one among other to test formulas for Bernstein-Sato poly- 
nomials of non-generic arrangements. A formula for the Bernstein-Sato polynomial of a 
generic hyperplane arrangement was given by Walther in 



Since the set E 



3.2. Topologically equivalent singularities 

Let f,g be two topologically equivalent singularities and assume that is known. 

|g27ria I = 0} is a topological invariant of the singularity {/ = 0} 



/ 



at the origin [TTIIIH] and every root belongs to (— n, 0) (Theorem 3.2), one can find an upper 
bound for bg{s) from the roots of &/(s) and use our algorithms for computing bg{s). The upper 
bound is constructed as n/3G£;('^ ~ Z^)' where E = {a + k\ aE Ej, k E Z,a + k ^ {—n, 0)}. 

In general it is complicated to check, whether two singularities are equivalent. However, 
there are some special families for which this can be done. This is the case of quasi-ordinary 
singularities, see e.g. [I6j. Let us see an example of a non-isolated one. 

Example 3.5. Let f = + x^y^ and g = f + x^y'^z. Since the corresponding discriminants 
with respect to z are normal crossing divisors, the associated germs at the origin define quasi- 
ordinary singularities. Moreover the characteristic exponents are in both cases the same and 
hence they are topologically equivalent, see e.g. P^ . 

The Bernstein-Sato polynomial of / at the origin has 27 roots, all of them with multiplicity 
one except for a = —1 which has multiplicity two. Here is the list in positive format. 



1, 



5 9 4 


13 2 3 


19 


5 11 17 


17 


11 


7 


19 


13 27 


7 21 9 13 5 


3 


7 


31 


7 


23 


6' 10' 3' 


10' 3' 4' 


20' 


12' 10' 12' 


20' 


12' 


10' 


12 


' 20' 20' 


6' 20' 20' 12' 4' 


2' 


12 


20' 


4 


' 20 
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The exponential of the above set has 24 elements. Each of them gives three candidates 
for bgfl{—s) except for —a = 1 which gives just two. For instance —a = 1/2 gives the 
following three possible roots. 

1 fl 3 5-1 

2 ^ l2'2'2i 

There are 71 possible roots in total. Note that using this approach we do not have any 
information about the multiplicities. Finally one obtains the roots for bgfi{—s). 

5 9 4 13 2 3 19 5 11 17 17 11 7 13 27 7 21 9 13 5 1 7 11 23 
' 6' 10'3' 10' 3'4' 20' 12' 10' 12'20' 12' 10' 20 ' 20 ' 6 ' 20 ' 20 ' 12 ' 4' 2 ' 12 ' 20 ' 20 



Observe that the Bernstein polynomials are very similar. The roots of bffi{—s) marked 
with a box have disappeared in bgfi{—s) and the ones in bold 3/2, 31/20 have become 
1/2, 11/20. In the table we put the information on timings in minutes: seconds format. Here, 
bf ct (computing Bernstein-Sato polynomial) and Sannf s (computing Ann£)[s](/*)) stand for 
corresponding procedures from the library dmod.lib. Also there are the minimal and the 
maximal time, spent for checking single roots. 



bfct(/) 


Sannf s((7) 


check all roots 


-a = 27/20 


-a = 1/2 


bfct(^) 


0:13 


0:45 


2:37 


0:04 


0:02 


33:03 



3. 3. A 'Campo 's formula 

The Jordan form of the local Picard-Lefschetz monodromy of superisolated surface sin- 
gularities was calculated by Artal-Bartolo in [2]. The main step in this computation was 
to present explicitly an embedded resolution for this family and study the mixed Hodge 
structure of the Milnor fibration. 

Since every root of the Bernstein-Sato polynomial belongs to the interval {—n, 0) (The- 
orem 3.2) and the characteristic polynomial is a topological invariant, using the results by 
Malgrange [T71 HB] , one can eventually provide an upper bound for the 6-function. Let us 
see an example that was not feasible even with the powerful specialized implementation by 
Schulze [3T] . 

Example 3.6. Let V be the superisolated singularity defined hj f = + {x'^z + + xy^). 
The characteristic polynomial is 

^ (t^-l)(t^-l)(t^^°-l) 

^ ' (t-l)(t30_l)(t24_l) • 

This polynomial has 76 different roots and thus we know in advance that the Bernstein- 
Sato polynomial (resp. the reduced one) has at least 77 (resp. 76) different roots. Using the 
above results in 230 possible candidates. Only 77 of them are roots of the 6-function indeed, 
all of them with multiplicity one. The total running time was 1 hour and 32 minutes. 

■1 27 101 41 17 83 103 43 53 29 107 23 89 109 71 91 37 73 31 

' 40' 120' 60' 20' 120' 120' 60' 60' 40' 120' 40' 120' 120' 120' 120' 40' 120' 40' 
113 37 47 19 77 97 39 13 49 59 79 33 119 3 4 121 47 161 181 61 
120' 60' 60' 20' 120' 120' 40 ' 20 ' 60 ' 60 ' 120 ' 40 ' 120 ' 5 ' 5 ' 120' 40' 120' 120' 60' 
71 27 91 41 143 163 61 21 73 83 31 127 49 167 187 43 149 169 131 
60' 20' 60 40' 120' 120' 40' 20' 60' 60' 20' 120' 40' 120' 120' 40' 120' 120' 120' 
151 57 133 51 173 67 77 29 137 157 59 23 79 89 139 53 179 6 7 
120' 40' 120' 40' 120' 60' 60' 20' 120' 120' 40' 20' 60' 60' 120' 40 ' 120 ' 5 ' 5 
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Remark 3.7. Spectral numbers are defined using the semi-simple part of the action of the 
monodromy on the mixed Hodge structure on the cohomology of the Milnor fiber [33], [36] . 
In [TTl Th. 3.3], |29l Th. 0.7] it is proved, that some roots of the Bernstein-Sato polynomial 
of a germ with an isolated critical point at the origin, can be obtained from the knowledge of 
the spectral numbers of the germ. Since spectral numbers do not change under /i-constant 
deformations, this also gives a set of common roots of the Bernstein-Sato polynomials, asso- 
ciated with the members of a /^-constant deformation of a germ. Therefore, they provide a 
lower bound for 6/(s), as well as an upper bound. 



4. Integral Roots of b-functions 

For several applications only integral roots of the 6-function are needed, e.g. [28]. We 
present here problems related to the so-called Logarithmic Comparison Theorem and In- 
tersection Homology D-module. Depending on the context local or global version of our 
algorithm is used. 

4-1- Upper bounds from different ideals 

Consider a left ideal / C Ann^[,] f. Then I + (f) C AnnD[s] f + (/) ^ D[s], that is 
the former is a proper ideal. Then define the relative 6-polynomial G K[s] to be the 
monic generator of (/+ (/)) nK[s], then | bj{s). Note, that quite often bf{s) = 0. But 
if 7^ 0, it gives us an upper bound for bf{s). In particular, one can take /, giving rise 
to a holonomic D[s]-module, that is GK.dimD[s]/I = GK.dimD[s]/ Ann^jjs] f^ = n + l. 

Since (s + 1) | | one can consider the reduced relative 6-polynomial 

bj{s) e K[s] to be the monic generator of (/+(/, ... , ^)) n K[s]. 

A prominent example of / as above is the logarithmic annihilator. Let / = Ann^j^j(/*) 
be the ideal in D[s] generated by the operators -P(s) G Ann/jf^] of total degree at most one in 



di. Let us define bf{s) := = {Ann^j^l^{f') + D[s]{f) ) nK[s]. The reduced bf^' 



is useful as well. 

4-2. Minimal integral root ofbf{s) and the logarithmic comparison problem 

Since every root of bf{s) belongs to the real interval (— n, 0), integral roots are bounded 
and therefore the whole Bernstein-Sato polynomial is not needed. Let us see an example 
that could not be treated before with the classical methods. 

Example 4.1. Let A be the matrix given by 

Xi X2 Xs X4 

A = \ x^ x-j Xs 

Xg Xio Xu X12 

Let us denote by Aj, i = 1,2,3,4, the determinant of the minor resulting from deleting the 
i-th column of A, and consider / = A1A2A3A4. The polynomial / defines a non-isolated 
hypersurface in C^^. Following Theorem 3.2, the set of all possible integral roots of bf{—s) 
is {11,10,9,8,7,6,5,4,3,2,1}. 
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Using the algorithm checkRoot with the logarithmic annihilator, see Section 4.1 above 



instead of the classical one, we have proved that the only integral root of b^p{s) in (—12, 0) is 
— 1. Hence —1 is the minimal integral root of bf{s). The following is the timing information 
of the whole procedure. Of course, —1 is always the root, but it is interesting to compare 
the timings of confirming this fact. 



Possible integral roots 


1 


2 ... 11 


Root ofb\^\s) ? 


Yes 


No 


Running time 


3:01:12 


~ 7:50 



This example was suggested by F. Castro- Jimenez and J.-M. Ucha for testing the Loga- 
rithmic Comparison Theorem, see e. g. |34j . 

The use of logarithmic annihilator allowed to reduce the computation time. However, for 
/ from this example it is known, that Ann£)^[s](/*) = Ann^^j_,j(/^) and this fact together with 
some homogeneous properties were used to compute other roots of bf{s), see Example 
below. 



4.4 



4-2.1. Quasi-homogeneous polynomials 

Assume F G Rn is a w-quasi-homogeneous polynomial with Wi ^ 0, that is, there are 
numbers Wi, . . . ,Wn such that with ^ = Yl^=i'^i^A has F = ^{F). Take c G K* and 
let us denote / = F\x^=i for some fixed k. We are interested in studying the relationship 
between the Bernstein-Sato polynomials of / and F. The result has been obtained working 
out directly with the functional equation. 

Proposition 4.2. Let F G Rn be a quasi-homogeneous polynomial with respect to the weight 
vector w = {wi, . . . ,Wn) . Assume Wk ^ for some k G {l,...,n} and define f to be 
the polynomial, resulting from making the substitution = c E K.* in F. Then bf{s) 
divides bpis). 

Proof. Consider the V^-filtration on Dn given by the variable Xk- Let P{s) G Dn[s] a dif- 
ferential operator satisfying the functional equation for F. There exists d > such that 
xf.P{s) G J2i>o^k ' ^0- From the homogeneity of F one can deduce that 

xA • F'+' = ^ + 1 - V WixA) • F'^+\ 

Let D' be the (n — l)-th Weyl algebra in the variables xi, . . . ,Xk, ■ ■ ■ ,Xn- Thus Vq = 
D'[xkdk] and xfP{s) ■ F*"*"^ can be written in the form Q{s) ■ F'"^^^ where the operator dk 
does not appear in Q{s) G -D„[s]. The functional equation for F has been converted in the 
following one: 

xtP{s) • = Q{s) • = xt bpis) • F' 
Now the substitution = c G K* can be made and the claim follows. □ 
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Example 4.3. The Bernstein-Sato polynomials of F = x'^z + and / = F\z=i = x"^ + are 



, / N / ( 5\ / 7\ ^ 4\ / 5 



6/ V 6/ V 3/ V 3 



hj{s) 

From the result of Kashiwara [12] one can see, blowing up the origin of F, that the last two 
factors are related to the 6-function of {z^ = 0}. This is a general fact. 

Example 4.4. Now, we continue Example |4.1 Let g be the polynomial, resulting from / by 



substituting xi, X2, x^, x^, x^, Xg with 1. Using Proposition 4.2 several times, one can easily 



see that bg{s) divides bf{s). Finally, the checkRoot algorithm is used to obtain that 

{s + l)\s + l/2)(s + 3/2)(s + 3/4)(s + 5/4) 
is a factor of bg{s) and therefore a factor of 

4-3. Intersection homology V -module 

In this part of the paper we introduce some new notation. We refer to [35j for further 
details. Let X be a complex analytic manifold of dimension n > 2, Ox the sheaf of holomor- 
phic function on X and Vx the sheaf of differential operators with holomorphic coefficients. 
At a point x G X, we identify the stalks Ox,x with the ring O = C{xi, . . . , x„} of converging 
power series and Vx^x with V = 0{di, . . . , On)- 

Given a closed subspace F C X of pure codimension p > 1, we denote by ifj^j(Cx) 
the sheaf of local algebraic cohomology with support in Y. Let X) C H^Y]i.^x) be 
the intersection homology Px-Module of Brylinski-Kashiwara. This is the smallest Vx- 
submodule of H^y]{^x) which coincides with ifp^j(Ox) at the generic points of Y . 

A natural problem is to characterize the subspaces Y such that C{Y, X) coincides with 
Indeed, from the Riemann-Hilbert correspondence of Kashiwara-Mebkhout, the 
regular holonomic Dx-module H^y]{'^x) corresponds to the perverse sheaf Cy[p], while 
X) corresponds to the intersection complex ICy- This way, the condition X) = 
H^Y]{^x) is equivalent to the following one: the real link of y at a point x G y is a ratio- 
nal homology sphere. Torrelli proved, that the following connection to local Bernstein-Sato 
polynomial exists. 

Theorem 4.5 (Theorem 1.2 in [35] )• -^^^ Y <Z X he a hypersurface and h G Ox,x o local 
equation ofY at a point y eY. The following conditions are equivalent: 

1. C{Y,X)y coincides with H^Y]{^x)y 

2. The reduced local Bernstein- Sato polynomial of h has no integral root. 

The proof of the theorem is based on a natural generalization of a classical result due to 
Kashiwara which links the roots of the 6-function to some generators of a G C. 
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z' + 



Example 4.6. Let Y be the affine variety in X = defined by the polynomial / 
{x'^z + y^){x'^ + y^z). The surface Y has the origin as its only singular point and thus the local 
6-function and and the global one coincide. The only possible integral roots are —2 and — 1. 



Now consider Jj, the Jacobian ideal of /, cf. Remark 2.8 Since the reduced Bernstein-Sato 
polynomial is required, the ideal 

AnnBw(r) + /^[s](/,J/,s + «) 



is used for checking rational roots, compare with Corollary 2.3 (2). We see that the above 
ideal is not the whole ring for a = 1 and hence the set of points x eY such that CiY, X)x = 
Hfy^iOx). isF\{0}. 

Using the implementation by Schulze [31] (based on Gauss-Manin connection), the com- 
putation of the whole Bernstein-Sato polynomial took 1236 seconds. While with our approach 
only 7 seconds were needed. 

Remark 4.7. Given Y as above, the set of points x G y for which the condition X)^ = 
H^Y]i^x)x is satisfied, defines an open set in Y that can be effectively computed with the 
stratification associated with the integral roots of the reduced local 6-functions, see the 
sequence of varieties ^ below. For instance, in Example 5.3, the open set is V{f) \ V{y, z). 



5. Stratification Associated with Local b-functions 



From Theorem 2.10, one can find a stratification of C" so that fe/,p(s) is constant on each 
stratum. The first method for computing such stratification was suggested by Oaku [26] (see 
also [25], [27] and [3] for further information). However, this method relies on the primary 
decomposition of commutative ideals. Following the ideas started in Section[2.2.2[ we propose 



a new natural algorithm for computing such a stratification. At first, a stratification for each 
root of the global 6-function is computed. Then one obtains a stratification, associated with 



the local 6-function, notably without any primary ideal decomposition, see Example 5.2 and 



5.3 below. We have created an experimental implementation, which was used for presented 
examples. The substitution of primary decomposition with elementary operations clearly 
decreases the total complexity of this algorithm. 



This is a natural generalization of Proposition 2.17 



Theorem 5.1. Let {Pi(s), 
and consider the ideals /„ 



0, 



1. Then one has 



,Pk{s),f} be a system of generators of AianDis]W) + D[s]{f) 
= (/ : (s + ay) + D[s]{s + a), for a root of bf{s) and i = 



rriaip) > i 



PG v(/«,nc[x]). 



Proof. Repeat the same argument as in Corollary 2.6| and proceed as in the proof of Propo- 
sition |2.17| using Lemmas 2.14 and 2.15 when necessary. □ 



Using the notation of the above theorem, let Va,i be the affine variety corresponding to 
the ideal Ia,i fl C[x]. Then 



a,mci — 1 



C---CK.nCK.,_i:=C", 



(6) 
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and rriaip) = iif and only if p G Va,i-i\Va,i- We call this sequence the stratification associated 
with the root a. Let us see two example^ to show how this result can be use to compute a 
stratification associated with local 6-functions. 



x^z^ -9/80y^z^ e C[x,y,z]. The global 



Example 5.2. Consider / = 
6-function is 

bfis 

Take Vi = V[x^ + 9/V - 1, 1^2 = Vix, y, z^ - 1) and I/3 = V{l9x^ + 1, Illy' - 80, z). 
Then V2 (resp. V3) consists of two (resp. four) different points and V3 C Vi, Vi fl V3 = 0. 
The singular locus of / is union of Vi and V2. The stratification associated with each root of 
hf[s) is given by 



[s + lf{s + 4/3)(s + 5/3)(s + 2/3). 

v2 



a 
a 
a 
a 



-1, 
-4/3, 
-5/3, 
-2/3, 



C V1UV2 
C V2UV3 

c V, c C3. 



c C3 
c C3 



From this, one can easily find a stratification of into constructible sets such that &/,p(s) 
is constant on each stratum. 



1 peC'\Vif), 

s + l peV{f)\{V,UV2), 

(s + l)2(s + 4/3)(s + 2/3) peV^\Vs, 

{s + lf{s + 4/3)(s + 5/3)(s + 2/3) p E V3, 

(s + l)(s + 4/3)(s + 5/3) peV2. 



Example 5.3. Let M be the surface in given by the polynomial / 
y^z — x^z'^. The global 6-function is 



x'^yz + xy"^ + y^ + 



bf{s) = (s + 7/6)2(s + 5/6)2(s + l)\s + 4/3)(s + 5/3)(s + 3/2). 

The singular locus of / is the union of V{x,y) and V{y,z). Now, consider the algebraic set 
V = Vi^x"^ + 4a;, xz, y, + z) = {0, gi, ^2}- Then the stratification associated with each root 
of bf{s) is given by 



a 
a 
a 
a 



-7/6, -5/6, 
-1, 

-4/3, -5/3, 
-3/2, 



C C3 
C C3 



C V{x,y,z) C V{x,y) C 

C V{y,z)cV{f) 

C V{x,y,z) = {0} 

C V = {0,qi,q2} C C^. 

Finally, we obtain a stratification associated with local 6-functions and the corresponding 
univariate polynomials. 



The hypersurfaces for examples have been taken from http://www.freigeist.cc/gallery.h.tinl 
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'l peC'\V{f), 

s + l peV{f)\{V{x,y)UV{y,z)), 

{s + + 7/6){s + 5/6) p e V{x, y) \ {0, gi}, 

bfM = \{s + lf veV{y,z)\{Q,q^}, 
(s + l)(s + 7/6)(s + 5/6)(s + 3/2) p = gi, 

(s + 1)2(5 + 3/2) p = g2, 

6/(s) p = 0. 

Remark 5.4. Note that one can define a stratification associated with the roots of the local h- 
functions, that is taking no multiphcities into account. We have observed that our algorithm 
is especially useful and very fast for computing this stratification. In particular, this is the 
case when each root has multiplicity one. Finally, also observe that in any case the global 
6-function is not actually needed, if an upper bound (or just a set containing the roots of 
bf{s)) is used instead. 

Remark 5.5. We see some common properties between the factorization of a Bernstein-Sato 
polynomial with the so-called central character decomposition by Levandovskyy [13j. In 
particular, for 6/(s) = nQeA('^ ~ a)™", where A C Q is the set of roots of &/(s), there is an 
algorithm for computing the following direct sum decomposition of the module 

D[s]/iAnno[s]in + if)) = ^M/(Ann^w(r) + (/)) : J(a)-, 

where J (a) = {bf{s)/{s — a)"^°'). We plane to investigate this topic further and provide cyclic 
D[s]-modules, corresponding to different strata. 

There is a very recent paper [23j by Nishiyama and Noro, where the authors build a strat- 
ification without using primary decomposition. The authors use initial ideals with respect to 
weight vectors in computations, which is a classical (cf. [28] alternative to the methods, utiliz- 
ing annihilators Ann£,[s]{f^). In pQ there is a comparison of performance of both approaches 
for the computation of Bernstein-Sato polynomials. Notably, no method is clearly superior 
over another. Rather there are classes of examples, where the difference is very distinct. 
In particular, initial-based method score better results on hyperplane arrangements, while 
annihilator-based methods are better at complicated singularities, which are not hyperplane 
arrangements. A comparison of two methods for stratification is very interesting and it is an 
important task for the future. However, it seems to us that the method we presented will 
allow more thorough analysis of the algebraic situation due to the applicability of central 
character decomposition. At the moment it is not clear, whether such a decomposition exists 
for initial ideals. 
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6. Other applications 



6.1. Bernstein- Sato Polynomials for Varieties 

Let / = (/i, . . . , /r) be an r-tuple in K[a;]^. Denote by K(S') the universal enveloping 
algebra [/(gt^), generated by the set of variables S = (s^), i, j = 1, . . . ,r subject to relations: 

[sij,Ski] = SjkSii — 6iiSkj. 

Then, we denote by Dn{S) := D„ (g)K K(5'). Consider a free lK[x, s, j]-module of rank one 
generated by the formal symbol and denote it by M = K[x, su, . . . , ' /^^ where 

= /f" ■ . . . ■ ff."\ The module M has a natural structure of left D„(S')-module. Denote 
by AnnD„{s){r) the left ideal of aU elements P{S) G Dr,{S) such that P{S) • = 0, that is 
the annihilator oi in Dn{S). 

Theorem 6.1 (Budur, Mustata, Saito [7j). For every r-tuple f = (/i,...,/r) € K[xY 
there exists a non-zero polynomial in one variable b{s) G K[s] and r differential operators 
Pi (5), ... , PriS) e Dn{S) such that 

r 

Y,Pk{S)fk-r = h{su + --- + Srr)-r e M. (7) 
k=l 

The Bernstein- Sato polynomial bf{s) of f = [fi, ... , fr) is defined to be the monic poly- 
nomial of the lowest degree in the variable s satisfying the equation ([T]). It can be verified 
that bf{s) is independent of the choice of a system of generators of (/i, . . . , fr). Then the 
Bernstein-Sato polynomials of / can be computed as follows 

(AnnB„(s>(/') + (/i, • • • , /r)) n K[sn + ■■■ + s„.] = (6/(sn + . . . + s„0)- 

In [1] an algorithm to find a system of generators of AnnD(^s){f^) was given. Moreover, in 
computing the intersection of an ideal with the univariate subalgebra an optimized algorithm 
(which uses linear algebra approach) was used. 



The above formula together with Theorem 2A_ can be used to check rational roots of 
Bernstein-Sato polynomials also for affine algebraic varieties. Hence, following Corollary 2.6 , 
a stratification associated with the local 6-functions can be computed. 

6.2. A remark in Narvdez's paper 

In [22], Narvaez introduces a polynomial denoted by f3{s) verifying Ann£)[s](/^) C 

Ann^j^j(/*). For all the examples treated in [22], he was able to compute an operator 

P'{s) G D[s] such that bf{s) — P'{s)f G Ann^j^j(/*). The last example in the paper is quite 
involved and could not be computed by using any computer algebra system directly. An 
iterated process for finding approximations of involutive bases was used instead. 
Indeed, for this propose the operator is not really needed, since 

bf{s)-P{s)feAnn%P^) ^ bf\s) = bf{s) ^ bf\s)\bf{s), 
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and thus after computing b^^\s), one only has to check whether each root of the latter 
polynomial is indeed a root of the 6-function and the same with the multiplicities. 
By definition the following inclusions hold 

m{^^ms]in + (/) ) C Ann«,j(/^) + (/) C Anno[s]{n + (/)• 

This implies that bf{s) \ b^f \s) \ I3{s)bf{s). Additionally, if /3(s) divides bf{s), then the poly- 



nomials b^j\s) and bf{s) both have the same roots and the above condition is equivalent to 



bf\s) and bf{s) 
ma{b^j\s)) = ma{bf{s)) for every root a of 

Example 6.2. Let / = (a;iX3 + a;2)(a;]' — Xg) be the last example from [22]. The Bernstein-Sato 
polynomial and the polynomial are respectively 

bf{s) = (s + lf{s + 3/4)(s + 3/8)(s + 9/8)(s + l/4)(s + 7/8)(s + l/2)(s + 5/8), 
= (s + 3/4)(s + 5/8)(s + l/2)(s + 3/8)(s + 1/4). 

Now one only has to check that all roots of /3(s) have multiplicity 1 as a root of b^}\s). 



This can be done using Theorem 



2.1 



with / = Ann^r^, (/"*) + (/). Using this approach the 



computations are trivial, less than 5 seconds. 



7. Conclusion and Further Work 

As we have demonstrated, the family of checkRoot algorithms (implemented in the library 
dmod. lib of SINGULAR) has many useful applications in the realm of Z)-modules. Nowadays, 
it is the only method that allows one to obtain some roots of the 6- function without computing 
the whole Bernstein-Sato polynomial. The latter is often infeasible despite all the recent 
progress in computational D-module theory. 

We emphasize, that presented techniques are elementary (by utilizing the principal ideal 
domain of the center K[s] of -Dn[s]) but very powerful from computational point of view. 
Many intractable examples and conjectures could be treated with this new method, as we 
have partially illustrated. Moreover, a stratification associated with the local 6-functions 
can be obtained without primary decomposition [26] as in the very recent [23]. It is very 
interesting to study these algorithms further and compare our approach with the one of [23] . 

Unfortunately, these techniques cannot be generalized for Bernstein-Sato ideals, since 
such ideals lie in K[si, . . . , s^] for m > 2. 

We have demonstrated that one can use the idea of checkRoot for checking rational 
roots of 6-function of a holonomic ideal with respect to a weight vector [28]. This gives an 
easier method for computing, among other, integral roots of such 6-functions, if an upper 
bound is known in advance. In this context, it would be very interesting to have a version of 
Kashiwara's result for some holonomic ideals and certain weights, since many algorithms in 
Z)-modules theory are based on integrations and restrictions which need minimal/maximal 
roots. 
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